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Abstract. We consider the problem of computing the family of operator norms recently 
introduced in [I]. We develop a family of semidefinite programs that can be used to ex- 
actly compute them in small dimensions and bound them in general. Some theoretical 
consequences follow from the duality theory of semidefinite programming, including a new 
constructive proof that there are non-positive partial transpose Werner states that are 
r-undistillable for arbitrary r. Several examples are considered via a MATLAB imple- 
mentation of the semidefinite program, including the case of Werner states and randomly 
generated states via the Bures measure, and approximate distributions of the norms are 
provided. We extend these norms to arbitrary convex mapping cones and explore their 
implications with positive partial transpose states. 



1. Introduction 

In [Ij we initiated the study of a family of operator norms that quantify the different 
degrees of entanglement in quantum states. Taken together with previous work [21 |3l 0], 
where special cases of the norms were considered, these norms have found applications 
to central problems in quantum information theory. Most importantly, to the problem of 
determining fc- entanglement witnesses and /c-positivity of linear maps, and to the existence 
problem for non-positive partial transpose (NPPT) bound entangled states. The work [1] 
in particular makes it clear that computing or bounding these norms even in special cases 
would have a significant impact on these problems. 

The primary goal of this paper is to continue the investigation of [1] . Here we focus on the 
development of algorithmic techniques to calculate and bound the operator norms, and we 
present further applications of the norms in quantum information. Specifically, we develop 
a family of semidefinite programs [SI El [7] that can be used to exactly compute them in small 
dimensions and bound them in general. Some theoretical consequences then follow from the 
duality theory of semidefinite programming, including a new constructive proof that there 
are NPPT Werner states [8] that are r-undistillable for arbitrary r [9]. We consider several 
examples via a MATLAB implementation of the semidefinite program. In particular, we 
show how they can be computed on Werner states and randomly generated states via the 
Bures measure [TUl [HI [12]) and provide approximate distributions of the norms. We also 
extend these norms to arbitrary convex mapping cones and explore their implications with 
positive partial transpose states [THl [HI [T3], and we apply them to a recent conjecture on 
the regularized relative entropy of entanglement |16| [T7] . 

The paper is arranged as follows. In Section [2] we present our notation and terminology 
and introduce the reader to the Schmidt rank of vectors and Schmidt number of density op- 
erators, /c-positivity of linear maps, and the Choi-Jamiolkowski isomorphism. The operator 
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norms are defined and some of tfieir most important properties are presented. In Section 2.3 



we will present an introduction of semidefinite programs and Section |3] will follow up by 
showing how they can be used to compute the operator norms in small dimensions and up- 
per bound them in general. We give a new constructive proof that there are NPPT Werner 
states that are r-undistillable for arbitrary r in Section |4] Section [5] contains MATLAB 
code that carries out the semidefinite programs, considers examples that demonstrate the 
performance of the semidefinite programs, and investigates the approximate distribution of 
the norms in small dimensions. The operator norms are then generalized to a larger family 
of convex cones in Section [6] and their properties are investigated. We conclude in Section [7] 
with further analysis of how the norms behave in the important case of projections. 



2. Preliminaries 

Our set-up is similar to that of [Ij and thus the interested reader is directed there to 
learn about concepts such as Hermicity-preserving maps, the Choi matrix of a map |181 ll9j. 
and the Schmidt number of density operators |20j. For us Tin denotes an n-dimensional 
complex Hilbert space and CiJ-Ln) denotes the set of linear operators acting on Hn- idn will 
represent the identity map on C{7in)- We will consider bipartite systems Tin Tim and 
assume that m < n. A unit vector \v) € 7i {a pure state) is denoted using Dirac bra-ket 
notation, with {v\ := \v)* . We will denote the computational basis vectors (i.e., the vectors 
with 1 in the i^^ component and in all other components) by 

If X G C{7i) is positive then we will write X > 0, and A > B indicates that A — B > 0. If 
it is important to note that a positive operator X is invertible, we will write X > 0. It will 
sometimes be convenient to denote the cone of positive operators in C{Ti.) by {C{7i))'^. A 
(mixed) quantum state is represented by a density operator p > that satisfies Tt{p) = 1. 
Whenever lowercase Greek letters like p or a are used, it is assumed that they are density 
operators. General operators will be represented by uppercase letters like X and Y. E := 
n Sij=o refers to the rank-1 projection onto the standard maximally entangled 

state. We will say that a Hermitian operator X = X* G C(7in) ® C{%m) is k-block positive 
(or equivalently a k- entanglement witness) if it is the Choi matrix {idn (E> ^)iE) of a linear 
map $ : C{'Hn) J~'{Ti-m) that is /c-positive. The scaled Choi matrix of the transpose 
map, S := {idn "X" T){nE), is a unitary that is referred to as the swap operator because 
5(|a) |5)) = |6) (g) |a) for any separable pure state \a) (8> \b). 

2.1. Relationship Between Schmidt Number and fc-Positivity. Observe that the 
set of states p with SN{p) < k is a closed convex cone if we remove the requirement that 
Tt{p) = 1. This cone will be denoted Sk- Given a convex cone C C C(7i), its dual cone is 
the convex cone defined through the Hilbert-Schmidt inner product as follows: 

C'^ := {X € C{n) : TviXY) > VF G C}. 

It is known that the dual cone of S^, the operators with Schmidt number no greater than 
k, is exactly the set of fc-block positive operators, and vice- versa |14l I15j . 

The following two theorems are each just a way of restating the fact that the cone of 
unnormalized states with Schmidt number at most k is dual to the cone of A;-block positive 
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operators. Theorem 2.1 in particular provides an important second characterization of k- 
positivity of a hnear map that is sometimes given as the definition of fc-positivity in quantum 
information theory 



Theorem 2.1. If ^ : ^ C{n n) is a linear map, then $ is k -positive if and only if 

{idn (g>'^){p)>0 \fpe C{rLn) «) C{rLn) with SN{p) < k. 
Theorem 2.2. If p & C{Hn) ® Cil-Ln) is a density operator, then SN[p) < k if and only if 
{idn ^ ^)(p) > Mk-positive $ : £('H„) >C('H„). 

These theorems are of theoretical interest, but are not of much practical use for testing 
fc-positivity or Schmidt number, since (for example) it is not possible to apply idn ® ^ 
to p and check positivity for every fc-positive map <1>. In some cases, however, an explicit 
finite set S of A;-positive maps is known for which [idn ® ^){p) ^ for all $ G 5 implies 
SR{p) < k. For example, if m = 2 and n = 2 or n = 3 then the transpose map T alone 
is enough to determine whether or not p is separable (i.e., SN{p) = 1) [21]. The fact that 
the transpose map can be used to determine separability in small dimensions has led to 
the study ol positive partial transpose (PPT) states [22|, which are density operators p such 
that {idn'^T){p) > 0. Throughout the rest of this paper, we will write the partial transpose 
operation {idn ® T){p) as p^ . 

The connections between /c-positivity and Schmidt number via dual cones have been 
studied quite a bit over the last few years. This basic theme of duality between /c-positivity 
and Schmidt number will be present throughout much of this paper. 

2.2. Family of Operator Norms. For 1 < k < m, the A;th operator norm introduced in 
[1] has the following form for positive operators X: 



X\\^^^^=snp[{v\X\v):SR{\v))<k} 
= sup {Tr(Xp) : SN{p) < A;}. 



The following two simple results about these norms were proved in [T]. Proposition 2.3 
shows that the problem of computing the norms for positive operators is equivalent to 
the problem determining /c-block positivity for arbitrary operators and is thus likely very 



difficult. Proposition 2.4 shows that we can nonetheless efficiently compute the operator 



norms when the operator under consideration has rank 1. 

Proposition 2.3. Let X G {C{'Hn) C,iT~im))^ be positive and let c G M. Then cl — X is 
k-block positive if and only if c> ||^||_5(^)- 

Proposition 2.4. Let p = \v){v\ G C{Hn) <8> C{7im) be a pure state and let {ai} be the 
Schmidt coefficients of \v) ordered so that ai > Q2 > • • • > ctm ^ 0. Then 



k 

\p\\s{k) = • 



1=1 
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2.3. Semidefinite Programming. Here we introduce the reader to semidefinite program- 
ming (SP), which we will see provides a step in the direction of being able to compute the 
operator norms defined above. Our introduction will be brief - for a more in-depth intro- 
duction and discussion, the reader is encouraged to read any of a number of other sources 
including [2'6\ |23 I27j . Most importantly, there are explicit methods that are able 

to approximately solve semidefinite programs of the type presented in this paper to any 
desired accuracy in polynomial time |5]. 

For our purposes, assume we have a Hermicity-preserving linear map ^ : C{'Hn) — >■ 
'C('Wm); two operators A £ C{7in) and B G C{T-Lm), and a convex cone C C [Cil-Ln))^ ■ 
Then the corresponding semidefinite program is given by the following pair of optimization 
problems: 

Primal problem Dual problem 

maximize: Th:{AX) minimize: Ti{BY) 

subject to: ^X) < B subject to: ^>t(y) > A 

X Y 

Though the semidefinite program ([T]) differs from the standard form of semidefinite pro- 
grams, it is equivalent and better suited to our particular needs. This form has been used 
very recently to solve other problems in quantum information [6l [7] . The interested reader 
is pointed to Appendix I for a discussion of how to convert between the form ([T]) and the 
standard form, as well as MATLAB code that performs the conversion in order to allow 
pre-existing software to solve these semidefinite programs of the form ([T| . 

We define the primal feasible set A and dual feasible set B to be 

A := {X e C : ^{X) < B} and B := {V e : {¥) > A} . 
The optimal values associated with the primal and dual problems are defined to be 
a := sup {TiiAX)} and /3 := mf^ {Tr{BX)}, 

and if ^ or ,S is empty then we set a = —oo or /? = oo, respectively. 

Semidefinite programming has a strong theory of duality. The theory of weak duality 
tells us it is always the case that a < (3. Equality is actually attained for many semidefinite 
programs of interest though, as the following theorem shows. 

Theorem 2.5 (Strong duality). The following two implications hold for every semidefinite 
program of the form ([T]) . 

1. Strict primal feasibility: If /S is finite and there exists an operator X in the interior of C 
such that ^{X) < B, then a = (3 and there exists Y £ B such that Tr(YB) = f3. 

2. Strict dual feasibility: If a is finite and there exists an operator Y in the interior of 
such that <I>"f(y) > A, then a = f3 and there exists X £ A such that Tt{XA) = a. 



There are other conditions that imply strong duality, but the conditions of Theorem 2.5 
(which are known as Slater-type conditions) will be sufficient for our needs. 

3. Bounding The Operator Norms 

Proposition |2.4| shows that we can compute the fcth operator norms of rank-1 opera- 
tors efficiently, since the Schmidt coefficients of a vector can be computed in O(n^) time. 
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However, Proposition 2.3 shows that the problem of computing these operator norms for ar- 
bitrary positive operators is equivalent to the problem of determining A;-block positivity for 
arbitrary operators and is thus likely very difficult. Here we develop a family of semidefinite 
programs that can be used to provide upper bounds on the norms in general and compute 
them exactly in low-dimensional cases. Additionally, some simple theoretical results that 
further establish the link between the kth. operator norm and fc-block positive operators 
will follow from the duality theory of semidefinite programming. 

Given a positive operator X £ (g) C{'Hm))'^ and a natural number k, we now 

present a family of semidefinite programs with the following properties: 

• Strong duality holds for each semidefinite program. 

• The optimal value a of each SP is an upper bound of \\X\\ ^^^^ . 

• There is an SP in the family such that the optimal value satisfies a = \\X\ 



\s(ky 



Let X G {C{'Hn)(E>C{'Hm))~^ be a positive operator for which we wish to compute \\X\ 
Let $fc : C(7im) C{T-lm) be a fixed /c-positive linear map and consider the following 
semidefinite program: 



Primal problem 

maximize: Ti{Xp) 



(2) subject to: 



mmimize: 
subject to: 



Dual problem 

A 

A/„ (g) /, 

y > 



> {idn<S)^l){Y) + X 



Tr(p) = 1 
p>0 

It may not be immediately obvious that this semidefinite program is in the form of 
([T]), so we first check that these problems are indeed duals of each other and form a valid 
semidefinite program. To this end, consider the linear map ^ : C{'Hn) ® C,{T~im) {i^{T~in)® 
C{n^)) e clUi) defined by 

" -{idn®^k){p) 
Tr(p) 



Then the dual map : (£(?^„) C{nm)) ® C{ni) £(^„) ® is given by 



Y * 
* X 



Finally, setting 



A:= X and B := 




1 



and recalling that the convex cone of positive semidefinite operators is its own dual cone 
gives the semidefinite program ^ in standard form. 

We now show that this program satisfies the Slater-type conditions for strong duality 
given by Theorem 2.5, It is clear that both a and /3 are finite, as Ti{Xp) < 



|X||andA>0. 

Both feasible sets are also non-empty (for example, one could take p to be any separable 
state, y = 0, and A > ||^||)- Strong dual feasibility then follows by choosing any y > 
and a sufficiently large A. Strong primal feasibility is not necessarily satisfied, however, as 
there is no guarantee that does not introduce singularities in p (for example, consider 
the zero map, which is /c-positive) . We could restrict the family of fc-positive maps that we 
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are interested in if we really desired strong primal feasibility, but strict dual feasibility is 
enough for our purposes. 

It follows from Theorem 2.1 that, for any fe-positive map the optimal value of the 
semidefinite program ([2]) is an upper bound of ||^||_5(-;j,-) ~ the supremum in the primal 
problem is just being taken over a set that is larger than the set of operators p with 
SN[p) < k. This leads to the following theorem. 

Theorem 3.1. Let X G (£(?^„) (g) CCHm))^ ■ Then 

\\'^\\s{k) ~ 111"'^ ~'~ ^11 • ^ ^"^ k-block positive^. 

Proof. Because is A;-positive if and only if is fc-positive, the dual problem Q can be 
rephrased as asking for the infimum of ||X + y || , where the infimum is taken over a subset of 
the A;-block positive operators Y G C{Hn) ® CiJ-Lm)- The preceding paragraph then showed 
us that 

ll-'^ll^^^-) ^ iiif {||^ + : Y is A:-block positive}. 

To see that equality is attained, choose Y := \\X\\^^^^I — X, which we know from Proposi- 
tion 2.3 is fc-block positive. Then 



|x + y|| - ||x + - ll^ll^^^j. 



□ 



In fact, it is not difficult to see that there is a particular fc-positive map such that 
is attained as the optimal value of the semidefinite program ^ corresponding 

to ~ simply let be the map associated with the operator — X via the 

Choi-Jamiolkowski isomorphism. 



3.1 



is that ||^||_5.(^) < + Y\\ for all 



One additional obvious implication of Theorem 

X G (CCHn) O C{Um)y and all fe-block positiveT^G CiUn) CiUm)- The following 
corollary shows that this can be strengthened into another characterization of A;-positivity. 

Corollary 3.2. Let Y = Y* £ C{T-Ln) (E> C{'Hm)- Then Y is k-block positive if and only if 
\\X\\sik) ^ 11^ + ^11 e (^(^") ® 



Proof. The "only if" direction of the proof follows immediately from Theorem 3.1 To see 
the "if" direction, assume that Y is not fc-block positive and choose X = cL — Y, where 
c G M is large enough that cl — Y > 0. Then, because Y is not fc-block positive, there exists 
a vector \v) with SR{\v)) < k such that (r;|y|f) < 0. Thus 

> {v\{cl - Y)\v) = c- {v\Y\v) > c = ||X + Y\\. 

□ 

Recall that if m = 2 and n = 2 or n = 3 then the transpose map T alone is enough to 
determine whether or not p is separable (i.e., SN(p) = 1 if and only if (idn tX" T){p) > 0). 
It follows that the semidefinite program ^ with k = 1 and $i = T can be used to 
compute for positive operators X. That is, the infinite family of semidefinite 

programs reduces to just a single semidefinite program in this situation. We can then use 
Proposition 2.3 to determine 1-block positivity of operators X G C{T-L^) ® C{7i2). 
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4. Undistillable Werner States 

Given a bipartite state /O, a natural question to ask is whether of not it can be transformed 
via local operations and classical communication to a maximally entangled pure state. If 
it can, p is said to be 1-distillable. It may happen that p itself cannot be transformed into 
a maximally entangled pure state in this way, but r copies of p can. In this situation, p 
is said to be r-distillable. If there exists an r such that p is r-distillable, then p is simply 
called distillable and otherwise p is said to be undistillable. The key connection between 
(un)distillability and the family of norms considered here comes from a result of [281 129] 
that says a state p is r-undistillable if and only if (/j®^)'" is 2-block positive. We will thus 
focus particularly on the 2-norm in this section. 

Not surprisingly, separable states are undistillable. It is also known that all states with 
positive partial transpose are undistillable [30' '3T] , but the converse remains an open prob- 
lem. That is, are there states with non-positive partial transpose (NPPT) that are undis- 
tillable? In this section we compute the 1-norm on a family of projections, and use that 
result to find states that are r-undistillable for arbitrarily large r. Note that this does not 
answer the question on whether or not NPPT undistillable states exist, however, because in 
our construction, the viable values of r are bounded above by a function of the dimension 
(n) of the Hilbert space. 

We shall focus on Werner states [8], which are a particular family of bipartite states in 
C{T-Ln) ® C{'Hn) that are central in quantum information. They are exactly the states that 
are invariant under any operator of the form U ®U, where U G C{T-Ln) is unitary, and they 
take the following form: 

Pa ■= —, — OiS) G C(T-Ln) ® /^(T-Ln) for some a G [—1, 1], 

n[n — a) 

where S is the swap operator as defined earlier. Werner states have become a subject of 
great interest because it has been shown that NPPT undistillable states exist if and only if 
there is an NPPT undistillable Werner state [32j. The following is a well-known result on 
Werner states. 

Lemma 4.1. If Pa' is r-undistillable, then pa is r-undistillable for all a < a' . 

In [T| it was shown that the Werner state p2/n (which has non-positive partial transpose) 
is r-undistillable if and only if the 2-norm of the projection onto the negative eigenspace 
of (p^^)^ is less than or equal to 1/2. More specifically, the projections of interest can be 
defined recursively as follows: 

nPl ■.= EeC{nn)®C{Hn), 

nPr-=[I-E)®nPr^l+E®{I-nPr-l) Vr>2, 

where / is the identity operator of appropriate size. Although the 2-norm of these projec- 
tions is still unknown, we can compute the 1-norm of each of these projections using the 
semidefinite program ([2]). 

Lemma 4.2. Let ^Pr be a projection as defined by the recurrence relations Q. Then 

1 1/ 2\'- 
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Proof. To see the ">" inequality, consider the separable vector \v) := |0) |0) G "H®^ 
^®2r-2 rpj^g^ ^gf^^g ^j^g quantity 

Cn,r := {v\nPr\v). 

It follows that 

c„,, = (0|(I - ^)|0)(0|„P,_i|0) + (0|S|0)(0|(/ - „P,_i)|0) 
n — 1 1 

— Cn,r—1 ~l~ (1 Cnr—lj 

n n 
2\ 1 

I Cn.r— 1 ~l~ 

nJ n 

1 1 / 1 2 



Standard methods for solving recurrence relations yields Cn,r = 2 ~ 2 ( ~ n ) ■ Noting that 



I "^'''15(1) — "^"'^ gives the desired inequality. 



To see the "<" inequality, we will use the dual form of the semidefinite program ([2]) with 
the transpose map ^i{X) := X'^ . To this end, notice that if A^^^' is the maximal eigenvalue 



of nPr 1 then (A™"^/ — nPr ) is positive semidefinite and so Theorem 3.1 says that 



II pII .^11 P J- ( \ 'TT-O'^ T pr \ r 1 1 1 1 \ maa; r II \max 

IP '"lls(l) — 1" \^n,r n-f^r / || ~ ||^7i,r || ~ ^n,r • 

In order to compute AJ^"^, let us consider the partial transpose of the family of projec- 
tions ([3]): 

^p[ = -sec{nn)(S)C{nn), 

n 

nP^ = -S®{I- nPr-l) + " -S) ® n^f-l Vr > 2. 

n n 
It is clear that the eigenvectors of nPr ^-re each of the form \x) ^ \y) for some eigenvector 
\x) of S and some eigenvector \y) of n-P^_i- If we recall that the eigenvalues of S are ±1, it 
follows that 

A-r = max{(l - -)A™-i + -, (1 + -jA^^i - -|. 

If K^r-i < \ then (1 + DK'^L^ - ^ < (1 - l)KT-^ + ^ < 2, so it follows via induction 
(and 'the fact that A^^f = ^ < 5) that A;^^^ = (1 - f )A^^^i + ^- We already saw that this 

recurrence relation has the closed form AJ^"^ — \~\\^~~] 1 which finishes the proof. □ 



2 2 \^ 

We can now state and prove the main result of this section. We use the above results 
together with results from [Ij in the proof. 

Theorem 4.3. Letn,r £N be such that p := ^}^'^^}_2y — ^/^ odd and a < 1(^+1) 

or if r is even and a < ^{ ''~\fp + 1); then the Werner state pa G C{'Hn) ^ C{T-Ln) is 
r-undistillable. 



Proof As ||n^'r||5(2) — 2 1 1 n Pr \\s{i) (see Theorem 4.13 of P), Lemma [42] above implies that 
(4) |UP,||5(2)<l-(l-^) . 
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The eigenvalues of (Pa^)^ are 

(1 — an)*" for m = 0, 1, . . . , r. 

In particular, (p^^)^ is nonsingular and has some positive eigenvalues, so the first half 
of condition (2) (with /c = 1) of Theorem 5.1 from ^ is satisfied. Now assume that 
P '■— ri^-(n-2Y — ^ that a = ^( ^''''''^l-^/p + 1). If we can show that pa is r-undistillable 



-(n-2y 

2 



then we are done by Lemma 4.1 Well, p > 1 implies that a > -, so the minimal positive 
eigenvalue A^j„ of (Pa^)^ is 1, and its maximal (in absolute value) negative eigenvalue X^ax 
is (1 -an)2r^/2l-i. We have 



n n\ y — (n — 2) 

Rearranging this expression yields, in the second equality. 



A;t... = i = («n-i)2r./2i-M(-^)'-i)>A^ 



iPr\\si2) 



2 1 ~ I — ' "tnax -1 II n II 

J I 1 - ||n-f>>||s(2) 



and where the final inequality comes from (|4]). Now by condition (2) of Theorem 5.1 from 
[T], we have that {p®^')^ is 2-block positive, and hence the result follows. □ 



Note that the value p of Theorem 



4.3 



is such that p > 1 if and only if n > ^ . 
Thus, for any r > 1, there is always some non-PPT Werner state that is r-undistillable 
as long as the dimension n is large enough. In fact, the dimension grows roughly linearly: 
is asymptotic to + 1- Also, if p > 1 then the result immediately implies that 

the a = 2/n Werner state is r-undistillable. Additionally, it is not difficult to see that if 
Pa S C{Hn)®C,{T~in) IS r-uudistillable then pa S C{T-im)®C{T~im) must also be r-undistillable 
for any m < n. Putting these facts together gives us the following slightly weaker (but much 
simpler) corollary of Theorem |4.3[ 



Corollary 4.4. If a < min{2/n, ln(2)/(r + 31n(2) - 1)} then pa S CCHn) CCHn) is 
r-undistillable. 

Similar results about r-undistillability of Werner states have appeared in the literature 
in the past. In [j9j it was shown that, for any fixed n > 3, there exist NPPT Werner states 
that are r-undistillable, though the region that was shown to be r-undistillable shrinks 
exponentially with r. Our result is stronger in that our regions shown to be r-undistillable 
shrink only like 1/r. On the other hand, for each fixed n our result only gives a region of 
NPPT r-undistillability for r < ln(2)(n - 3) + I. 

5. MATLAB Implementation and Examples 

The semidefinite programming method of bounding the operator norms has been im- 
plemented in MATLAB, with the script available for download from the website http:^ 



//www. nathanieljohnston. com/ index. php/schmidtoperatornorm/. In order to test the 



semidefinite programs, we will need a theoretical result to compare the computed results 
to. To this end, we compute analytically the fcth operator norms of Werner states. We also 
look at the operator norms of randomly generated states from the Bures measure. 
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5.1. Werner States. Here we compute the kth operator norms of Werner states. Specifi- 
cally, the following result shows that if a < then each norm coincides with Loa . If q > 



then 



Po 



IS(1) 



is smaller, but the rest of the norms are all equal to \\po 



Proposition 5.1. Let pa G C{T-Ln) 'S> C{'Hn) be a Werner state. Then 
II II l + |min{a,0}| , n n 1 + la 



\Pa\ 



5(1) 



and 



Po 



\S(k) 



n{n 



a) 



for 2 < k < n. 



n{n — a) 

Proof. Throughout the proof, we will work with the operator '■= n{n — a)po 
to simplify the algebra. To see the result when a < 0, note that for any k, 

W^aW^fi.-, = \\l - oiS\\g^f^^ < ||/|L,,^ - allS'lL,,,^ = 1 



aS 



a, 



\S{k) II" "'^llS'(fc) - II IIS(A:) "iniS(fc) 

where the inequality comes from the triangle inequality and the rightmost equality comes 
from the fact that = 1, which is easily verified. To see the other inequality, choose 

|f) := |0) (8) |0) and observe that 

n-l 

{v\X\v) = m ^ (0|)(/ - aS){\0) |0)) = l-aJ2 (0|i)(i|0)(OU)(i|0) = 1 - a. 

i,j=0 

On the other hand, if a > 0, then for any vector \v) = \a) <E> \b) , it follows that 

{v\X^\v) = {{a\ {b\){I - aS){\a) = 1 - a{{a\ \a)) = 1 - a\{a\b)\'^ < 1. 

Furthermore, equality can easily be seen to be attained when \v) = |0) (E) |1), which shows 
that ||Xq,||^j.^^ = 1. To see the result for /c > 2 and a > 0, use the triangle inequality again to 

< 1+Q. To show that equality is attained, let Iv) = ^(|0)(g)|l)-|l)(g)|0)) 



and observe that {v\X, 



1 + a. Since \v) has SR{\v)) = 2, the result follows. 



□ 



The performance of the semidefinite programs for the 1-norm is analyzed in Table [TJ If 
the transpose map is used, then we know that the semidefinite program must give exactly 
1 1 '^"11 5(1) ^^^^ n = 2, which it indeed does. In fact, the map 4>i defined by ^i{X) = 
Tt{X)I — X that is used as the basis of the reduction criterion [32] also gives the correct 
answer in this case. For n = 3, the transpose map still happens to give the correct answer, 
though the reduction criterion map gives a strict upper bound when a > 0. 









Upper bound computed using... 


n 


a 


Exact Pa 


$i(X) = X^' 


$i(X) = Tv{X)I-X 


2 


1/2 


1/3 


0.3333 


0.3333 


2 


-1/2 


3/10 


0.3000 


0.3000 


3 


1/2 


2/15 


0.1333 


0.2000 


3 


-1/2 


1/7 


0.1429 


0.1429 



Table 1 . The exact 1-operator norm of various Werner states as well as the 
computed upper bounds obtained by using the semidefinite program defined 
by one of two different positive linear maps. 
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5.2. Randomly Generated States Via The Bures Measure. As another example, 
we consider random density operators distributed according to tlie Bures measure [TOl [H] , 
which can be generated quickly via the method of |12j . We then investigate the general 
behaviour of the kth operator norms of a density operator in C{'H2) ® ^^2) and C{'}i^) ® 
C{T-L^) relative to its eigenvalues. 

In particular, Figure [T] shows how the 1-norm is distributed compared to the two largest 
eigenvalues A3 < A4 in £('^2) ® C{T-L2), based on 2 x 10^ randomly-generated density op- 
erators. It is not surprising that the 1-norm lies between the A3 and A4, since A4 is equal 
to the 2-norm, and it was shown in [1] that the {n — l)-norm in C(T-Ln) (S" C(T-Ln) is always 
at least as big as the second-largest eigenvalue. We see that the 1-norm typically is much 
closer to A4 than A3. 




Figure 1. Approximate distributions of the 1-norm and the two largest 
eigenvalues of random Bures density operators in £(^2) (X" £(7^2) • 

The 1-norm in this case was computed using the semidefinite programming method of 
Section [sj A similar plot was presented in [33j for what was called the "maximum local 
eigenvalue" , which coincides with the 1-norm for positive operators. There it was similarly 
observed that the 1-norm typically lies closer to A4 than A3 under the Hilbert-Schmidt 
measure. 

Figure [2] shows how the 1 and 2-norms typically compare to the two largest eigenvalues 
As < Ag in C{T-L^) (>S> C{T-L^), based on 10^ randomly-generated density operators. As before, 
it is not surprising that the 2-norm lies between As and Ag. However, it was shown in fT] that 
there exist density operators p G J~-{7i3) ® C{Hz) for which A5 < ||/o||5(i) < Ae. Nonetheless, 
this situation seems to be extremely rare, as ||/o||5(i) generally lies between As and Ag. 
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Value 



Figure 2. Approximate distributions of the 1 and 2-norms, as well as the 
two largest eigenvalues of random Bures density operators in C{T-L^) (E) C{'H^) . 

Because the semidefinite programming method of Section [3] does not produce the exact 
value for the 1 and 2-norms in C{'H^) C{'Hs), the values of the norms used for Figure [2] 
are estimates that were derived from a simple genetic algorithm. 

6. Norms Restricted to Other Convex Cones of Operators 

We will now see that many of the results for the kth. operator norms actually hold in the 
much more general setting of arbitrary convex mapping cones of operators. We will begin 
by defining the notion of a mapping cone, which was originally introduced by Stormer|13J. 

Definition 6.1. Let S C C(C{'Hn) , ^{Tim)) be a cone of completely positive linear maps. 
S is said to be a mapping cone if ^ o ^ ^ S and ^ o ^ ^ S whenever $ G 5 and ^ is 
completely positive. 

Mapping cones appeared recently in [15] as a way of generalizing the dual cone relation- 
ships between A:-block positive operators and operators with Schmidt number no greater 
than k. These dual relationships can be seen implicitly in the semidefinite programming 
results of Section [3| so it is no surprise that the notion of mapping cones provides a natural 
generalization in this setting as well. Mapping cones can be defined without the restriction 
that they be a subset of the completely positive maps, though the definition provided will 
be better-suited to our purposes. 

We will say that a cone of operators C C (£(^„) (g) C{T-lm))~^ is a mapping cone if the 
cone of associated linear maps (via the Choi-Jamiolkowski isomorphism) is a mapping cone. 
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If necessary, we will specify whether we mean a mapping cone of operators or a mapping 
cone of linear maps, but our meaning should be clear from context. 

Definition 6.2. Let X G {C{%))^ he positive and let C C [C{'Hn) ® C{%rn))^ be a closed 
convex cone. Then we define the C-operator norm of X , denoted ||-^||^, by 

\\X\L :=sup|Tr(X/))). 
^ pec J 

It is easy to see that this defines a valid norm if C is a mapping cone. It is also a norm for 
many other convex cones of interest - all that needs to be checked is that C contains a full 
set of n'^m? linearly independent operators. Observe also that if C = 5^ and X > then 
this definition reduces to exactly The norm H-^^Hj^ has a similar interpretation to 

that of the kth operator norms as well. We can think of ||^||^ as roughly measuring how 
close X is to an operator in C. 

It is trivial to see that if C Q "D, where P is another closed convex cone, then < 
In particular this implies that ll-'^H^ < ll-'^H always. Additionally, several of the 
characterizations of the fcth operator norms carry over in an obvious way to this more 
general setting. 

Proposition 6.3. Let X G {C{%n) ® ^(^m))^ be positive. Then cL — X if and only 

ifc>\\X\\^. 

Proof. By definition, c/ — X G if and only if 

TY((c/ - X)p) = c - Tv{Xp) > 'i pee. 
This if true if and only if c > ||X||^, completing the proof. □ 

Now let X G {C{T-in) ® ^iT~{-m))^ be positive and consider the following semidefinite 
program. 

Primal problem Dual problem 

, , maximize: Tr{Xp) minimize: A 

^ ' subject to: Tr(p) = 1 subject to: A/„, >Y + X 

p€C Y eC^ 

It is easy to see that these problems are indeed duals of each other and form a valid 
semidefinite program, using the same method as was used in Section [3] to show that the 
semidefinite program ([2]) is valid. Strong dual duality also holds in this setting. The main 
difference here is that we have p G C and Y G rather than p, y > - we could have 
stated the semidefinite program ^ in terms of the cone 5^, but then it would become less 
clear how to actually implement the semidefinite programs and compute upper bounds of 
ll"'^lls'(fc) ^s™S ^-positive maps. 

Just as in the case for the A;th operator norms, the theory of semidefinite programming 
leads to the following two results. We state them without proof, as their proofs are almost 



identical to the proofs of Theorem |3 . 1 1 and Corollary 3.2, respectively. 
Theorem 6.4. Let X G (£(^n) ® C{nm))^ be positive. Then 

||x||^ = inf {||x + y|| :yGC^}. 
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Corollary 6.5. Let Y = Y* e CCHn) ® >C(^„). Then Y e if and only if 

\\x\\^ < \\x + Y\\ yx e {CiUn) ® cCHm))^. 

6.1. Application to PPT States. Given any positive linear map <I> : C{l-Ln) C{'Hm), 
there exists a natural convex cone C$ associated with 

C$ := {X G iCCHn) CiUm))^ : (idn > O}. 

Given any such convex cone, || • ||^^ is indeed a norm and we are able to compute 
to any desired accuracy via semidefinite programming, as seen in the previous section. In 
fact, ll^ll^ is exactly what is computed by the semidefinite program It follows that 



k 



\\'^\\s{k) ~ ^^^*fe mI^IIcj. ' /c-positive} . 

In the case of the transpose map T : C{Tin) CiTin), Ct is exactly the cone of unnor- 
malized PPT states, and so the norm || • ||^^ can be seen as a measure of how close a given 



6.4 



operator is to having positive partial transpose. It is known [M] that the dual cone of the 
PPT states is given by 

= ^X = X* e CCHn) CCHm) :X = Y + ZioT some Y>0, {idn ® T){Z) > o}. 

This leads immediately to the following characterizations of ||p||^^ via Theorem 

Proposition 6.6. Let p £ C{'Hn) C{T-Lm) he a density operator. Then 

\\p\\^^=\ni{\\p + Y\\:{idn®T){Y)>Q]. 

7. Norms on General Projections and a Conjecture of Brandao 

We have seen that the A;th operator norms of orthogonal projections have several ap- 
plications within quantum information theory. One more reason for studying these norms 
comes from their appearance in a conjecture of Brandao P^, which asks whether or not 
there exists a < e < 1 such that, for all n and all orthogonal projections P = P* = G 

C{nn)®C{'Hn), 



/ N II „ii /rank(P) 

(6) II^Ld) > V^^- 

In order to answer this question, recall the inequalities (6) and (7) of Theorem 4.13 of p^. 

For any fixed n G N and rank(P) < n, then Inequality (6) of implies that the statement 
of Inequality ^ holds for some < e„ < 1. Similarly, if rank(P) > n then Inequality (7) of 
[1] implies that the statement holds for some < 1- Thus, the statement of Inequality ^ 
holds in any fixed finite dimension. 

Nonetheless, Inequality Q can be seen not to hold as n tends to infinity when rank(P) ~ 
n via methods of convex geometry. In particular, we prove the following result using the 
ideas presented in [Tj 



Theorem 7.1. There exists a universal constant C, independent of n and k, such that for 
a general orthogonal projection P G C{'Hn) ^ C{T-Ln) with rank{P) < kn, we have 

h k 

-<\\PUk)<c-. 
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1-norm lower bounds 




Figure 3. A comparison of the various lower bounds for the norm on pro- 
jections in the n = m = 15 and k = 1 case. Taken together, Inequahties 
(6) and (7) from [1] show that Inequahty ([6| holds for all fixed, finite 
(Inequality ^ is shown with e = 0.99). 



n 



Proof. The left inequality is true for all projections simply by Theorem 4.13 of [T]. We will 
prove the right inequality by making use of the "tangible version" of Dvoretzky's theorem 
that appears in [IT]. First note that for any orthogonal projection P £ CiTin) ^ C^Hn), 



(7) 



\P\ 



S{k) 



sup \ . 2, '^i '■ {o^i} 3.re the Schmidt coefficients of \v) \. 

\v)eRange{P) ^ 



This characterization of ||-P||s(fc) appeared in the proof of Theorem 4.15 of [T] and also in [3], 
so we use it without proof. But now by associating Tin ^ Tin with C{'Hn), the quantity ^ 
equals 

(8) 



K 

sup |. ^s^(^) : \\A\\f = l,si{A) > S2{A) > ••• > s„(A) > singular values of a\, 
\ i=i 

where R is the subspace of C{l-Ln) associated with the range of P through the standard 
bipartite vector to operator isomorphism. So now the goal is to show that there exists 

a constant C such that y^X^iLi sl{A) < Cy^^V^||^||f for A in general subspaces R of 
dimension kn. To this end, we need to bound the constants h and M of Dvoretzky's 



theorem. It is trivial to see that \J^i=iS~{A) < \\A\\f and that equality if attained for 
some operators A, so 6 = 1. 

To upper-bound M, recall from ^17j that the expectation of the operator norm, E||^ 
upper-bounded by ^ for some absolute constant Cq. Thus 



is 



M :=E 



^sj{A) < E{VksiiA)) = VkE\\A\\ < Co 
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It follows via Dvoretzky's theorem that there is a constant c such that if we choose 
e = l/{Co^/c), then for general subspaces R with dim(i?) < ce^CQkn = kn, we have 



k 

. ^ s^{A) < (1 + e)M\\A\\p < (1 + ^ 



^11 



□ 



In the case when k = 1, Theorem 7.1 tells us that for general projections of rank n, 



||-P||5(i) < C/n, so Inequality ([g]) can not hold as n tends to infinity if e < 1. However, In- 
equality ([G]) is still relevant as it only needs to hold for certain projections to have important 
implications. If it did, it would imply that the regularized relative entropy of entanglement 
[6A\ I35| is super-additive [31|. This in turn would imply that QMA{k) = QMA(2) for all 
k > 2 via a result of Aaronson et. al. [_37j. 
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Appendix I: Implementing Semidefinite Programs 

The presentation of a semidefinite program in Section [2^ was as an optimization problem 
defined by a Hermicity- preserving linear map <I> : C{'Hn) ^(^m), two operators A G 
C{T-Ln) and B G jO{Tim), and a closed convex cone C C C{T-Ln) with the following primal and 
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dual forms: 



(9) 



Primal problem 

maximize: Ti{AX) 
subject to: <^>{X) < B 
X £C 



Dual problem 

minimize: Ti{BY) 
subject to: ^\Y) > A 
Y eC^ 



Here we will show explicitly, in the special case of C = {C(7in))~^ , how to convert the above 
semidefinite program into the so-called standard form of a semidefinite program, defined by 
a vector c ^ Hi and operators D G i^{T~ip) and {Fi] G C{'Hp), with the following primal 
and dual forms: 



(10) 



Primal problem 

maximize: x*c 

subject to: Yli=i ^i^i < D 



Dual problem 

minimize: Tt{DY) 
subject to: Tr{FiY) = Ci yi<i<£ 
Y >0 



Once the conversion from form ^ to form (10) has been carried out, the problem can 
be given to a semidefinite program solver to be solved. In particular, we provide a MAT- 
LAB front-end that carries out the upcoming conversion and uses the SeDuMi semidefinite 
program solver |38| to compute the solution. 

Thus, assume that you have a semidefinite program in the form Define a linear map 



MX) 






-X 



Then the requirements that ^{X) < B and X > are equivalent to the single constraint 



B 




The dual map ^'^^ : {C{%rn) © C{'Hn)) — ^ C{%n) acts on block diagonal operators as 



Y 
Z 



Thus, the semidefinite program ^ can be written in the following form: 
(11) 



Primal problem 

maximize: Th:{AX) 
subject to: '^{X)<D 



Dual problem 

minimize: Tr(DVF) 
subject to: ^^{W) = A 



where D 



B 




and W 



Y 




Note in particular that we can replace the inequal- 
ity in the dual problem ([9| by equality in ( 11 ) because of the flexibility that was introduced 
by the arbitrary positive operator Z. Now let Ea and Fa be families of left and right general- 
ized Choi-Kraus operators for ^ (that is, operators such that "^{X) = Y^^EaXFa). Denote 
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the {g, /i)-entry of X by Xgh and the (i, j)-entry of Ea and Fa by Saij and faij, respectively. 
Then 



where 



^'(X) = EgXFa = ^ y^eaigXghfahjj.. = XghGgh 



a gh 



gh 



Ggh ■— Yl 



(ialgfahl 
^a2gfahl 



(^algfah2 
ea2gfah2 



^a(m+n)gfahl (^a(m+n )gfah2 



(^alg fah(m+n) 
^a2g f ah{m+n) 

^a{m+n)g fah{m+n) 



Then by examining the equahty constraint in the SDP (11), we see that, for all g,h, 



Yl,i fagiWilCalh Y^i fagiWi2ealh 

fagiWiiea2h Ylii fagiWi2ea2h 



X/i fagiWi(^m^n)^alh 
y^j j fagiWi(jfi^n)^a2h 



fagi'U^il^a{m+n)h Xyj fagiWi2G-a{rn+n)h 
'^AGhgW). 



(m+n)^a{m+n)h^ 



It follows that the semidefinite program (11) can be written in the form: 
(12) 

Primal problem Dual problem 

maximize: vec(X)*vec(A) minimize: Tr(Z)VF) 

subject to: Xl^^ XhgGhg < D subject to: Ti{GhgW) = agh V 1 < 5, /i < n 

1^ > 

where vec(X) and vec(j4) are the vectorizations of X and A, respectively, that are obtained 
by stacking the columns of the matrices on top of each other into a column vector in the 
usual way. The semidefinite program (12) is in standard form, so it can now be input into 
a semidefinite program solver. 

This transformation of a semidefinite program of the form ([T]) into a semidefinite pro- 
gram in standard form has been implemented in MATLAB as a front-end for the SDP 
solver SeDuMi. The code and usage instructions can be downloaded from |http : / / wwwT] 
[nathaniel j ohnston . com/ index . php/ quantumsedumi/ , 

'^Department of Mathematics & Statistics, University of Guelph, Guelph, ON, Canada NIG 
2W1 



^Institute for Quantum Computing, University of Waterloo, Waterloo, ON, Canada N2L 
3G1 



